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Abstract 

We study initial boundary value problems for linear scalar partial diflterential equa- 
tions with constant coefficients, with spatial derivatives of arbitrary order, posed on 
the domain {t > 0, < x < L}. We first show that by analysing the so-called global 
relation, which is an algebraic relation defined in the complex fc-plane coupling all 
boundary values of the solution, it is possible to identify how many boundary condi- 
tions must be prescribed at each end of the space interval in order for the problem 
to be well posed. We then show that the solution can be expressed as an integral in 
the complex fc-plane. This integral is defined in terms of an i-transform of the initial 
condition and a t-transform of the boundary conditions. 

For particular cases, such as the case of periodic boundary conditions, or the case 
of boundary value problems for second order PDEs, the integral can be rewritten as an 
infinite series. However, there exist initial boundary value problems for which the only 
representation is an integral which cannot be written as an infinite series. An example 
of such a problem is provided by the linearised version of the KdV equation. Thus, 
contrary to common belief, the solution of many linear initial boundary value problems 
on a finite interval cannot be expressed in terms of an infinite series. 

1 Introduction 

An evolution partial differential equation in one space dimension is characterised by its 
symbol, which we denote by Lo{k). This means that a particular solution of the equation is 
given by 

E{x,t) = eikx~^{k)t^ 
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Physically significant examples of scalar evolution equations are: 

(a) the Schrodinger equation with zero potential 

iqt + qxx=0, uj{k)=ik'^; (1.1) 

(b) the diffusion-convection equation 

qt^qxx + Qx, uj{k) ^ ~ ik; (1.2) 

(c) the Stokes equation 

qt + qxxx + qx = ^, uj{k) ^ i{k - k^). (1.3) 

We note that equations (|1.1() and (|1.3I) are the linearised versions of the nonlinear Schrodinger 
and of the Korteweg-deVries equations respectively. 

Notations 

(i) qo{x) denotes the given initial condition, and qoik) is the Fourier transform of qo{x): 

qQ{x) = q{x,0), < X < L- qo{k) = / qo{x)e~'^'^'' dx, keC. (1.4) 

(ii) {fj{t)}Q~^ and {gj{t)}o~^ denote the boundary values of the solution at a; = and 
X — L respectively, while {fjityk)}^^^ and {gj{t,k)}Q~^ denote certain t-transforms 
of fj and gj: 

f,{t)^diq{Q,t), t>0; fj{t,k)^ f e'^^'^'fj{s)ds, fc G C, t > 0, (1.5) 

Jo 

gj{t) ^ diq{L,t), t>0; gj{t,k) ^ [ e'^^'^^'^j (s)ds, keC.,t>0, (1.6) 

where j = 0, . . . , n — 1. 

(iii) The domain D in the complex fc-plane is defined by 

D = {keC: Rew(fc) < 0}. (1.7) 

D'^ and denote the part of D in the upper and lower half of the complex fc-plane 
respectively: 

D+ = {keC:ke D, Im(fc) > 0}, D- ^ {k e C : k e D, Im{k) < 0}. (1.8) 

The oriented boundary of £)+ and is denoted by 91)+ and dD^, where the ori- 
entation is such that the interior of the domain D is always on the left of the positive 
direction. 
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Statement of the problem and assumptions 

Let q{x,t) satisfy a linear evolution equation with symbol uj{k) in the domain 

{t>Q, 0<x<L}, 

where L is a finite positive constant. We assume that Lj{k) is a polynomial of degree n such 
that the equation uj{k) — has n distinct roots, and that Reuj{k) > for fc £ M. Assume 
that the initial condition qo{x) is a given, sufficiently smooth, function. 

(i) Determine the number of boundary conditions that must be prescribed at x = and x — L 
in order to define a well posed problem. 

(ii) Given appropriate boundary conditions at the two ends of the space interval, assuming 
that these given functions have sufficient smoothness and are compatible with go (x) at x = 
and X — L, construct the solution q{x,t). 

The classical approaches 

The form of the particular solution E{x, t) suggests that the most convenient representation 
is the one obtained by an x-transform. 

1. An x-transform 

Separation of variables gives rise to an 7i-th order differential operator, acting on the interval 
< x < L. 

• n = 2 

For second order differential operators there exists a well known spectral theory which states 
that if one boundary condition is prescribed at each end of the interval, then it is always 
possible to expand a function in terms of the complete set of the associated eigenfunctions 
Thus for PDEs of second order, using separation of variables it is always possible 
to solve boundary value problems on the finite interval using generalised Fourier series. 
For example, it is elementary to solve equation with Dirichlet or Neumann boundary 
conditions in terms of a sine or a cosine series respectively. Similarly, the Dirichlet problem 
for equation (|1.2(l can be solved in terms of a generalised Fourier series over the complex 
eigenvalues ^ + |, m £ Z. 

• n = 3 

For third order differential operators, there exists a complete basis of eigenfunctions only for 
particular types of boundary conditions These do not include the uncoupled boundary 
conditions arising in the typical boundary value problems for third order evolution PDEs. 
Thus it is not possible to solve equation with uncoupled boundary conditions by an 
x-transform. Furthermore, it is not possible to use the spectral theory of the associated 
differential operators Lf = /"' + /' to decide how many boundary conditions must be 
prescribed at each end of the interval in order to guarantee that the resulting boundary 
value problem is well posed. 

2. The Laplace transform 

The particular solution E{x, t) can be rewritten in the form e^*+*'^(*)^, where k satisfies the 
n-th order equation 

s + ijj{k) = 0. 
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• n = 2 

Evolution PDEs with spatial derivatives of up to second order can be solved by a Laplace 
transform. For example, for equation H1.2|) one can construct a representation involving the 
exponential 

gSt+i[l+(l+4s)i/^]f 

However, it is easier to construct the solution using the associated x-transform. 

• n = 3 

The application of the Laplace transform in this case fails. Consider for example equation 
First, it is not immediately clear how many boundary conditions must be imposed at 
each end of the interval to guarantee well posedness. Furthermore, even if one prescribes one 
boundary condition at a; = and two aX x = L (which is the correct number of conditions), 
one obtains a representation for the Laplace transform q[x, s) which has singularities along 
the real line, see Appendix (when n = 2, the singularities in the complex s-plane are along 
the half line arg(s) = (p, with 7r/2 < </? < 37r/2). 

The new method used here to analyse boundary value problems for equations with spatial 
derivatives of arbitrary order is mathematically straightforward, still it yields results which 
are either very difficult or cannot be obtained by the standard approaches. This method, 
which is the implementation to this class of problems of the general approach introduced by 
one of the authors 1", involves the steps outlined below. 

Reformulation of the PDE 

A given evolution PDE with symbol Lij(k) can be written in the form 

where q{x, t) is a solution of the PDE, and the function X(x, t, k) is given by the formula 

X{x,t,k) ^^Cj(k)diq{x,t), (1.10) 

where the coefficients Cj(fc) are known polynomials in k. For example, for equation 

n = 2, ci(fc) = —k and C2(fc) = i. The explicit form of Cj{k) for an arbitrary ujik) is given 

in section 121 

The global relation and well posedness 

Equations p.9|l - l|1.10|l . and Green's theorem applied to the domain {[0,L] x [0,i]}, imply 

n-l 

5^c,-(fc)(/,(t,A:)-e-'*^^5j(i,fc)) =9o(fc)-c"W*g(t,fc), fc G C, (1.11) 



where goi fj and gj are defined in the notations (equations (|1.4|l - H1.5|l land q{t,k) denotes 
the Fourier transform of q{x,t). 
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Equation is the so-called global relation. The crucial observation is that the functions 

fj and gj depend on k only through Lj{k) It follows that these functions are invari- 
ant under any transformation of the complex fc-plane that leaves iv{k) invariant. These 
transformations are determined by the roots of the equation 

uj{k) = uj{X) given hy Ao(A;) = fc, Ai, . . . , A„_i. (1.12) 

Replacing k by X{k) in equation (|1.11|) we obtain a system of n equations 

n-l 

^ c,{Xi{k)) (^f,{t, k) ~ e-'^'('=)^5,(t, k)) = q,{Xi{k)) - e-('=)*<j(t, Xi{k)), l = 0,..,n~ 1. 

3=0 

(1-13) 

Ignoring for the moment the unknown function q{t,k), equations (|1.13() can be considered 
as n equations coupling the 2n unknown functions {/j, 5j}o^^- Therefore we expect that if 
n boundary conditions are prescribed, hence that n of the functions {fj{t,k), gj{t,k)} are 
known, the remaining n unknown functions can be obtained by solving this system of n 
equations. However, not all choices of n boundary conditions yield a well posed problem. 
Indeed, it turns out that the contribution of the additional unknown function q{t, k) vanishes 
if and only if a set of N and n — N boundary conditions are prescribed at a; = and x = L, 
respectively, where N depends in general only on the symbol a;(fc): If n is even, then 
N = n/2, while if n is odd then N is either {n+ l)/2 or (n — l)/2 depending on the sign of 
the highest derivative. 

We will show in sectionOthat N can be computed as follows: the set of prescribed boundary 
conditions must be such that when solving for the unknown functions in the set {fj(t, fc)}, 
the expression involving q{t,k) must be bounded as A: — > oo when k € , and when 
solving for the unknown functions in the set {gj{t, k)}, the expression involving q{t, k) must 
be bounded as fc — > oo when fee . 

The representation of the solution 

(i) For problems where an x-transforms exists, equation Hl.ll|l porvides an alternative way 
of constructing the solution. Consider for example equation (|l.l|l and assume for simplicity 
that homogeneous Dirichlet boundary conditions are prescribed. Then /o = go ~ and 
Ai = —k, hence equations H1.13|l become 

if\^ie~^kL~^ = qo(k)^e"'"'q{t,k), 

^/i-^e^'^^Si = qo{-k)-e'''"'q{t,-k). 

Subtracting these equations we obtain 

Evaluating this equation at the values of k for which the coefficient of gi vanishes, i.e. 
k — kjn :— m £ Zi, and using the definition of (?o(fc) and q(t, k) we find 

g^^™^ {t, fc„0 = e-''=^*<7^^'") (fc™), (fc,„) = sm{kr,x)q{x)dx. (1.14) 

Jo 
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Equation ^.14^ ) can be inverted by the well known formula and then q(x, t) is expressed in 
the form of a sine series. 

Similarly, consider equation H1.2|l and assume for simplicity that homogeneous Dirichlet 
boundary conditions are prescribed. Then proceeding as in (i) above, in analogy with 
(|1.14|l . we find 

(7(i,fc,„)-g(i,i-fc,„) = 0-^'^'"-''='")* (go(fcm) - go(j - fcm)), ^™ = + rneZ. (1.15) 

The left hand side of equation 11.15h ) involves the difference [e~'*^™^ — e~'*^*~'^™)^] . Using 
the orthogonality of this expression with respect to e**^"^ — e'*^*"*^"-*^, n =i m, equation (|1.15b ) 
can be inverted, yielding 

li^.t) = ^ E - e'^*"'"^") c-C^"-*")* (go(fc„0 - qo{^ - km)) ■ (1.16) 

(ii) For problems where there does not exist an x-transform, the alternative approach also 
fails. Consider for example equation with g(0,i) = q{L,t) = qx{L,t) = 0. In this case, 
substituting in the global relation /q = 50 = ffi = 0, and evaluating the resulting equation 
at fc, Ai(fc) and A2(fc), we obtain a system of three algebraic equations for /i, /2 and §2. 
The determinant of this system vanishes at fc = km, m G Z, where for large fc, arg(fcm) is 
either 7r/6, 5n/6 or 37r/2 (see section 0)) The functions /i, /2 and 52 involve e*''''"''" ^'^j and 
for large k these exponentials become unbounded at fc = km- Hence, we cannot evaluate this 
system at the zeros of its determinant. 

It will be shown in section 01 that the solution of an arbitrary boundary value problem for 
an evolution equation with symbol Lij(k) can be expressed in the form 

„ ri— 1 

9(a:,i) = ^{l e^'=^--W*go(fc)& - / e^'=---W*^c,(fc)/,(t,fc)dfc 




9^+ j=o 



n-1 



^k{x — L)—ijj{k)t 



J2c,{k)~g{t,k)dk\ , (1.17) 



where qo, dD~^,dD~ are defined in the notations and {fj,gj}a^^ can be computed as follows. 
Let fp{t) and gr{t) denote the prescribed boundary conditions, where p and r take N and n — 
N integer values respectively. Then fp and gr can be computed immediately (see equations 
1)1. 5|l and 1)1. 6|l ). Let fp{t) and gR{t) denote the remaining unknown boundary values. 
Solving the n algebraic equation 1)1. 13(1 for fp and gp it follows that these unknown functions 
can be expressed in terms of go, {/p, ffr}, and of a term involving q{t, A(fc)) and 1/A(fc), where 
A(fc) denotes the determinant of the relevant system. Because of our characterisation of well 
posedness, the expressions of q{t, A(fc)) appearing in the representation of fp are bounded 
as fc ^ 00, fc G therefore, using Cauchy theorem in the domain it follows that 

these expressions give a zero contribution if A(fc) ^ 0, fc S D^. If A(fc) has zeros in 
then the contour dD'^ must be deformed to pass above all these zeros. The contribution of 
the zeros can alternatively be written in terms of the associated residues. The analysis of 
the terms gp, is similar. 
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Regarding equation we note that the associated determinant A(fc) has no zeros in D. 

Actually it will be shown in section^ that the solution of the boundary value problem for 
this equation with the boundary conditions 

g(0,t) = /oW, q{L,t)^go{t), q,{L,t) ^ gi{t), t > 0, (1.18) 

is given by 



1 f°° „ ^ f ■! 

q(x, t)^— e'^"-*('=^'= ^^oik)dk + — / e'kx-^(k-k _ k^)foit, k)dk 
Joo 27r Jqjj+ 

+ ± [ e'^i-^)-^i^-^')\l-eyg^{t,k)+ikUt.k) (1.19) 

+ ^ / e--(--^)*..A(,fc) + /.(,fc).fc + ^ / e-(^-)-(-^^)*^.(i,.).fc, 

where the contours dD^ and 91)^ are shown in figure ^ and the functions /i, /2 and g2 
are given in terms of the given initial and boundary conditions by 

*/i = ^ [e-''=^(7V(A2,i)-iV(Ai,<))+e-^^i^(iV(fc,i)-iV(A2,i)) + e-''^^^(iV(Ai,<)-iV(fc,t))] , 

/2 = ^ [e-''=^(A2Af(Ai,i)-Ai7V(A2,i))+e"''^^^(fciV(A2,t)-A2iV(fc,t)) 
+e-'^^^(Ai7V(fc,t) -fc7V(Ai,t))] , 

52 - -)-[(7V(fc,i)(A2-Ai) + 7V(Ai,i)(A2-fc) + 7V(A2,i))(fc-Ai)], 
A(fc) 

where A(A:) and N{k, t) are defined by 

A(fc) = e-''=^(Ai - A2) + c-'^^^[k - Ai) + c-'^i^(A2 - k), (1.20) 

N{k, i) = (1 - fc2)(/o(i, A:) - e'^^-^.go(t, A:)) - ifce-^'=^5i(t, fc) + 9o(fc), (1-21) 
and Ai, A2 are the two roots of equation (|1.12|) . 



Advantages of the method 

The method introduced by one of the authors 0j does not rely on separation of variables. 
Thus this method can be used to obtain appropriate integral representations of the solution 
of PDEs, for which the spectral theory of the associated linear differential operator fails. 
From the above discussion, it appears that the advantage of the new method presented here 
manifests itself mainly for equations of order greater than two. However this method has also 
certain advantages for the classical second order equations, particularly in the case of more 
complicated boundary conditions, such as Robin conditions. This is illustrated in section 
15.21 and further discussed in sectional 
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Figure 1: The domain D and the contours dD^ (in red) and dD (in green) for equation 
p.3|l . Note that has two connected components. 



2 A reformulation of evolution equations 

It can be verified (see 6 ) that a PDE with symbol uj{k) can be written in the form 

|^g-»fcx+a.(fe)i^^^ _ |^g-»fcx+c^(fc)t^^ =0, keC, (2.1) 

where the function X(x, t, k) is given by the formula 

uj(k) ~ Lu(~idr) , ^ , , 

X{x,t,k)=i ^ ^ > q{x,t). 2.2 

k + lOx 

2.1 Examples 
(a) The equation 

The symbol is 

ijj{k) — ik^ . 

Indeed, it can immediately be verified that e**-'^"*'^^* satisfies this equation. Hence 

X(x,t,k)=i—-——-q^ — — q = iqx~kq. (2.3) 

k + tOx k + lOx 

This PDE is thus equivalent to the expression 

(e-'^^+'^^'q^^- (e-''''=+'''^\iqx-kq)^ =0, fc G C. (2.4) 
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(b) The equation (fTT^ 

The symbol is 

ijj{k) = fc^ — ik. 

Therefore 

, ,^ . -ik + k'^ + dl + dx (k + idx)(k - idx - i) 

^' = * fcTTaf = k\idx — = ^ ^'^ + ^ ^ ^ 

This PDE is thus equivalent to the expression 

{e-'^'^+^^^-'^^'q^^- {e-'^'-'+^^^-'^^\{l+ik)q + qx)^ =0, fc G C. (2.6) 



(c) The equation ()1.3|l 

The symbol is 

uj{k) ^i{k-k^). 

Hence 

, . i(fc-fc3)-(9, + 9g) _ (fc + idx -k^- id'x) 

= + ^ fc^^g ->q=ik^ -l)q-ikqx-qxx. (2.7) 

This PDE is thus equivalent to the expression 

^^_,kx+^ik-k')t^^J^ _ |^^-^kx+^(k^k^)t^^p _ _ ^^^^ _ ^ ^ ^ (3.8) 

3 Well posed problems 
3.1 The global relation 

We first derive the global relation for the case of equations 
Example (a) 

Since X{x, t, k) = iqx{x, t) — kq{x, t), it follows that the global relation is 

zhit, k) - kfoit, k) - e-^^^ {lUt, k) - fcgo(i, k)) = qo{k) - e''='*(j(i, k). (3.1) 
Example (b) 

Since X{x, t, fc) = (1 + ik)q + qx, the global relation is 
(1 + tk)Ut, k) + f,{t, k) ~ e-^^-^ ((1 + ik)Ut, k) + hit, fc)) - go(fc) - e('='-*^-)*g(t, A:). (3.2) 
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Example (c) 

Since X{x, t, k) — (fc^ — \)q — ikq^ — Qxxi the global relation is 

{k^ - l)/o - ikh - h] - e-''^^ [(fc2 - - ikg, - 52] = qo{k) - S'^-'^'^'qit, k). (3.3) 

3.2 The analysis of the global relation 

The roots of the equation w(fc) = i^(A), which are defined in terms of k (see equation (|1.12() '). 
can be distinguished by fixing their asymptotic behaviour: 

A,„(fc) ~ e^'^™/", k^oo, l<TO<n-l. 

For the examples (a) - (c) the roots different from k are given below: 

(a) iX^ -ik"^ ^{X-k){X + k) = {), Xi = -k. 

(b) X^ -iX-k^ +ik= iX-k)iX + k -i) = 0, Xi^i-k. 



(c) (A-A3)-(fc-fc3) = (fc_A)(A2 + Afc + P-l)=0, 

A2 2 2- 



Before turning to the examples (a)-(c), we make the following observations, valid for an 
arbitrary Lu{k): 

(i) The functions qo{k) and q{t, k), which are entire functions of k, are bounded as fc — > 00 
in C~, the lower half of the complex fc-plane. 

(ii) The functions e^^^q^ik) and e''''^q{t, k) are bounded as fc — > cxd in C+, the upper half of 
the complex fc-plane. 

(iii) The functions fj and gj, which are entire functions of fc, are bounded as fc — > 00 if 
Rew(fc) < 0. 

Example (a) 

In this case, since 

Rew(fc) = Re(ifc2) ^ -2Re(fc)Im(fc), 

it follows that 

D = {fc e C+ : Re(fc)Im(fc) > 0}. (3.4) 

Hence D is the union of the first and third quadrant in the complex plane. Since Ai = — fc, 
equations (|1.13l) become 

ih{t, fc) - kfoit, fc) - e-*^-^ {tgi{t, fc) - k~go{t, k)) = qo{k) - e^'^'qit, fc), 
i/i(t,fc) + fc/o(i,fc)-e*'=^(i.gi(t,fc) + fc.go(i,fc)) = qo{-k) - e'^'^'qit, -k). (3.5) 

We consider equations H3.5|l as two equations relating the four terms /o, /i, .go, 3i, treating 
for the moment the function q{t, fc) as a known function. The dependence of the solution of 
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the system (|3.5|l on the terms involving q{t, k) is the crucial factor that determines its well 
posedness. 

To explain how to compute N , we distinguish two illustrative cases. 
(a.l) q(0,t) =fo(t), q(L,t) = go(t) 

In this case, the two functions /o(i, k) and go{t, k) are known, and we must solve the system 
H3.5|l for /i and gi. We concentrate only on the dependence of the solution on the term 
q{t, k). This yields 



ifi 



(known) + e*''^* 
(known) + e*''^* 



q(t,k)-q(t,-k) 



ikL 



l-c 
q(t,k) 



2ikL 



ikL 



q(t,^k) 



1 



o.2ikL 



The terms involving the function k) in equation Ij^.fih .) are hounded in k as k 

fee D+ = Dr\£+. 

Indeed, for k in C+, 



JlikL 



q(t,k) - q(t, ~k) 



JlikL 



1 



,2ifcL 



q(t, k) - q(t,-k), 



(3.6) 

CXD if 



and both these terms are bounded in C+. Since the term e is bounded in D, our claim 
follows. 

Similarly, the terms involving the function q(t, k) in equation Ijtj.Bb ) are bounded as A; — > oo 
if fc e D- =DnC-. Indeed, 



ikL 



q(t,k)~e"'^q(t,~k) 



1 - e 

and both terms are bounded in 



2ikL 



e-^'^q(t,k) 



^ — ikL 



q(t, -k), 



(a.2) q(0,t) =fo(t), qx(0, t) = fi(t) 

In this case we obtain a system for the two functions go(t,k) and gi(t,k). Considering 
explicitly only the terms involving the function q(t, fc), we find 



kgo 



known + e**^^* 



JlikL;, 



(known) 



ik-'t 



J2ikL 



q(t,k) +q(t,-k) 

2QikL 

q(t, k) + q(t,~k) 



2e 



IkL 



In this case not all the terms containing the function q(t, k) are bounded as fc — > oo for 
fc e D~ ^ Dr\C~. For example, 



r,2ikLA 



q(t, fc) + q(t, -fc) 



r^ikL 



q(t, fc) 



-ikL 



2e 



,ikL 



q(t, -fc) 



and the second term is bounded in C^, but the first term is bounded in C^. Thus in the 
above example, as expected, iV = 1, i.e. one boundary condition must be prescribed at each 
end of the interval. 
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Figure 2: The domain D and the contours dD^ (in red) and dD (in green) for equation 



Example (b) 

The domain D is in this case given by the interior of the hyperbola (see figure 

D^{keC: Re(fc)2 - lm{kf + Im(fc) < 0}. (3.7) 

D+ — D n C+ is the interior of the upper branch of this hyperbola, and = D Ci is 

the interior of its lower branch. 

Since Xi ^ i — k, equations yield 

(1 + ^fc)/o(^, k) + hit, k) - c-*^^ [(1 + ik)Ut, k) + gi{t, k)] - qo{k) - e^^' g(t, k), 

-tkfoit, k) + Mt, k) - e-^t'-'^)^ i-tkgoit, k) + .gi(i, k)) - qo{i ~ k) ~ e('='-^'=)*g(t, z - k). 

(3.8) 

As for the previous example, the analysis of the boundedness of the terms containing the 
function q(t, k) implies that one boundary condition must be prescribed at each end of the 
interval, hence N — 1. This analysis is similar to the analysis performed in the case of 
Example (a). 

Example (c) 

The domain D is (see figure^ 

D = {ke C : Im (k) [3Rc (k)'^ - Im (kf - l] < 0}. (3.9) 
Evaluating equations l|1.13|l at k , Ai and A2, we find 

{k^ - l)/o - ikh - h 



-e-'^'^lie ~l)~g,~ik~g,^~g2] - qo{k) ~ e'^^-^"^'q{t,k), 
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(A? - l)/o - zAi/i - /2I - e-'^i^ [(A? - 1)50 - ^1,91 - 32] = <7o(Ai) - e*('=-'=')*g(Ai, i), 



(A^ - l)/o - zA2/i - h 



-iX2L 



[(^2 - 1)50 - iA2gi - 02] 



'?o(A2)-e*("-"')*(7(A2,i). 

(3.10) 



It turns out that, for a well posed problem, one boundary condition must be prescribed at 
X = and two boundary conditions must be prescribed at x = hence A'^ = 1. 
We illustrate this by two particular examples. The first set of boundary conditions, with 
= 1, yields a well posed problems. The second, with N = 2, does not. 

(c.l) q(0,t) = fo(t), q(L,t) = go(t), qx(L,t) = gi(t) 

In this case the unknown functions /i, /2 and §2 are given by 

ifi = (known terms) 



+e 



-,-iAiL / 



(m - qiX2)) + e-»^^^(g(Ai) - q{k)) + e-^'^^qiX,) - g(A2)) 



A{k) 



f2 = (known terms) 



+e 



i{k-k^)t 



-iXxL ( 



(kq(\2) - \2q(k)) + e-'^^^(Aig(fc) - fcg(Ai)) + e-*'^-^(Aiq(A2) - A2^(Ai)) 



A(fc) 



ff2 = (known terms) 



+e 



i{k-k^)t 



q(Xi)(k - Ai) + g(A2)(A2 - fc) + g(fc)(A2 - Ai) 
A(fc) 



(3.11) 



where A (A:) is given by equation (|1.2U|) and to simplify the notation we have suppressed the 
t dependence in q(t, fc). 

The terms containing the function q(k) in equations H3.11b .b') are bounded m fc as fc — s- 00 
if fc £ = D n C"*". This follows from the observation that if fc G then Ai and A2 
are in C~. For example, consider the function fi(k). Since e*'''^, e^^^^^ and e"*"^^^ are all 
bounded for fc in C+, the bracket appearing on the right hand side of i/i, as fc 00, is 
asymptotically given by 

(g»(fc-Ai)L _ p^(fc-A.)L)^(fc) _^ _ c'ik~X2)L^(^Xi) + (e'(k-X^)L _ 1^(^X2) 

Ai — A2 

All terms in this expression are bounded when fc G C+. In addition, Q^(k-k^)t bounded 
for all k G D, and the claim follows. Similarly, the terms containing q(k) in equation H3.11b ) 
are bounded m fc as fc — > oo if fc £ D~ — D n C~. 

(c.2) q(0,t) - fo(t), q(L,t) = go(t), qx(0,t) = fi(t) 

In this case, the unknown functions /2, gi and 52 are given by 



/2 = (known terms) 



igi — (known terms) 



aikL I 



(A2 - Ai)g(fc) + e'^=^(Ai - k)q(\2) + e'^i^(fc - A2)g(Ai) 



A(fc) 



(3.12) 
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-i\iL I 



{q{k) - q{\2)) + e-^^^^(9(Ai) - q{k)) + e-^'=^(g(A2) - g(Ai)) 



A(fc) 



52 = {known term,s) 



i(k-k^)t 



-*^^^Ai(q(A2) - q{k)) + e-'^-^^Mm - q(Ai)) + C-'^^k{q{M) - qiM) 



A(fc) 



where A(fc) is given by 
A(fc) 



ikL I 



(A2 - Ai) + e*^^^(Ai ~k) + e*^i^(fc - A2) 



iXiL ( 



As for example (a. 2), not all the terms containing the unknown function q(t, k) in equations 
(|3.12|l are bounded for all k £ Z?+ or fc G . As an example, consider the terms in H3.12b ). 
which should be bounded as fc ^ 00 for all fc £ D~ . Choose k such that A(fc) £ and 
A2(fc) G D+. Then e"-^^ and e*^i^ are not bounded as fc ^ co, while e*'''^^ is bounded. Since 
asymptotically Ai ~ e^^^^'^fc, it is easy to verify that, for k G such that Ai(fc) G -D^, 
the dominant term in the denominator is e*'''^^. Hence (ll-{.12b ^. as fc ^ 00, is given by 

^.ik-k-)t [e-2^>ii(g(fc) „ g(A2)) + e*'=^(g(Ai) - ^(fc)) + e'^^^<?(A2) - <7(Ai)] 

and the term e'^^^q{k) ~ qI-Hl-x) bounded for fc G (to simplify the notation we 

have again suppressed the t dependence in qit, fc)). 

4 The integral representation of the solution 

We first derive equation p.l7|l . The global relation is 



<?(t,fc) = e--«*<7o(fc)-e 



o{k)t 



f{t,k)-e-'''^~g{t,k) 



where 



/>, fc) = ^ c,{k)f,{t, fc), g{t, fc) = ^ c,{k)g,{t, fc), 



q{t,k) 



-ikx 



q{x,t)dx, qQ{k) — q{0,k). 



Taking the inverse Fourier transform of q{t, fc), we obtain 



/ ^ 1 
Ztt 



ikx — i 



^W*{go(fc)-/>,fc)+e-^'=^.g(i,fc)} 



dk. 



It follows from the definition of D that for fc ^ Z) and for any t > 0, the functions 
^-uj{k)tj(j-^ fc) and e'^^^'^^g{t, fc) are bounded as fc -> 00. Thus 



Qikx~u>(k)tj(^l^-^ is analytic and bounded for fc G C+ \ Z?, 

Qik(x-L)-i^(k)t~(j^-^ is analytic and bounded for fc G \D. 
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An application of Cauchy theorem thus yields 



e 



dD+ 



3 JdD- 

and ((TT7|) follows. 

We now derive the solution representation for our three examples. In the case of the two 
second order examples, we will not need to use the general integral representation as a series 
representation can in these case be deduced in a more direct way. However, the integral 
representation can alternatively be derived, and we will show in the next section how it can 
be reduced to the one obtained here. 

4.1 Example (a) 

Consider equation with the boundary conditions 

g(0,t) = /o(i), q{L,t)=g^{t), t>0. (4.1) 

Equations (|3.5() become 

z/i - e-'^^^iffi = N{k,t)-e'''"'q{t,k), 
ifi-e'^'^tgi ^ N{-k,t)-e''"'q{t,-k), (4.2) 

where N(k,t) is the known function 

N{k, t) = k (^fo{t, k) - e-*'=^5o(i, fc)) + go(fc). (4.3) 
Subtracting equations H4.2|l we obtain 

^^^kL _ e-'kL^ig^ ^ ^^^^ _ _ e»fc'*(q(i, k) ~ q{t, ~k)), km G M. (4.4) 

We evaluate this equation at the values of k for which c'*^^ — q^^^l _ j 



TOTT 

L 



The definition of q{t, k) implies 



si\Y{kmx)q{x, t)dx = ^e [iV(fc™, t) - N{-~km,t)] 



Hence 

g(x,t) = ^^sin(fc„x)e-^'="*[7V(fc™,i)-^(-fc™,i)], k^^^- (4.5) 
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4.2 Example (b) 

Consider equation (|1.2|l with the boundary conditions (|4.1I) . Equations (|3.8|) become 
/i-e-''=^gi - iV(fc,t)-e('='-^'=)*g(i,fc), 

(4.6) 

where N(k,t) is the known function 

iV(fc, i) = -(1 + zk) (/o(fc) - e-^^^5o(fc)) + go(fc). (4.7) 
Subtracting equations H4.6|l we find 

(e-'^'-''^^ - e-'''^^ gi = N{k, t) - N{i ~ k,t) ~ e^k^-'''^* (git, k) - q{t, i - k)) . (4.8) 
We evaluate this equation at the values of k for which Q-^i^-f')^ — ^-ikL _ 

TTLTT I 

k = k^^ — + -, me Z. (4.9) 
1j z 

Equation 14.8|) and the definition of q{t^ k) yield 



(e-^fc™- _ e-'(*-fc™)-) q{x, t)dx = e-C^^-''^™)* (iV(fc„, t) - N{i - km, t)) . (4.10) 



This equation can be solved for q{x, t) using orthogonality, i.e. the fact that if fc„j is defined 
by then 



(^Q-ik^x _ ^-i{^-krr.)x\^ (^^^k„x _ gi(*-fc„)x^ ^ 2L5mn- (4.11) 

Thus equation (|4.10() yields 

q{x, = ^ E e-^'""^'^'")* (e^'^"" - e^(^~'=™)-) (Af(fc™, i) - N{i - fc™, t)) . (4.12) 

m 

4.3 Example (c) 

Consider equation (|1.3|l with the boundary conditions (|1.18() . Equations H3.10|l become 

- f2-ikfi+c-"''^92 = N{k,t)-e'(-''-''"'>'q{t,k), 
-/2-*Ai/i+e-*^i^52 = iV(t,Ai)-e'('=-'=')*g(t,Ai), 
-/2 - «A2/i + e-*^^^52 = iV(t,A2)-e'('=-'=')*g(t,A2), (4.13) 
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where the known function N{k,t) is given by p. 2111 . In this case, the representation (|1.17|l 
becomes 



dD+ 

M--L)-ik'-^'^)\(k^ - l)g, - ik~g^ - g,)dk \ , (4.14) 



q{x, t) = i- I / e*'=^-'('=-'=')*go(fc)dfc - / c''^---'('=-'=')*((A;2 - l)/o - zk'h ~ f2)dk 
e 

dD- 

where Z?+ is the interior of the upper branch of the hyperbola (|3.9() . and D~ is the exterior 
of its lower branch (see figure ^ . 

Solving equations (|4.13|) for /i, /2 and 52 and substituting the resulting expressions in 
equation 14.14|l we find equation (|1.19|) . We emphasise that the solution of the system 
H3.10|l includes also the terms appearing on the right hand side of equation H3.11|l . However 
these terms do not contribute to the solution. Indeed, after multiplication by cT'^i^^^ )* the 
terms in (|3.11b .b) are analytic and bounded as fc — > 00 in , while the terms in (|3.11b ) are 
analytic and bounded for fc — > 00 in D~ . Then an application of Jordan's lemma implies that 
these terms gives zero contribution provided A(fc) ^ for k in D. This is indeed the case, 
as it can be proven rigorously, using the fact that this determinant is an entire functions of 
finite order (the order of an entire function is a measure of its rate of growth as /c — > 00 - in 
this case the order is in fact equal to one). For such entire functions there exists a rich and 
classical theory, starting with the Hadamard factorization theorem, which implies that such 
a function has infinitely many zeros accumulating at infinity. The determination of these 
zeros in the complex plane is more difficult: in the present case, the general theory implies 
that the zeros must lie "near" the lines {arg(fc) — 7r/6, arg(fc) = 57r/6, arg(fc) — 37r/2}, at 
least asymptotically. By using the special symmetry properties of A(/c) one can then prove 
that asymptotically they lie exactly on these lines, and that they are always outside D. For 
the proof, which is tedious but elementary, see pii| . 

5 Alternative representations for second order prob- 
lems and more complicated boundary conditions 

5.1 Alternative representations 

The general representation (|1.17|) implies alternative representations for the Dirichlet bound- 
ary value problems associated with the equations (|l.l|l and H1.2|l . Regarding these alterna- 
tive representations, we note the following: 
(i) For example (b), equation (|1.17(l becomes 

q{x,t) = — (/ e''=^-('='-*'=)*go(fc)'^fc - / e'^'=-^'^^-''^^\{l + ik)fa + fi)dk 
27r iJ-oo JdD+ 

+ ^k)~g^ + ~g,)dk^ , (5.1) 



dD- 



where and D are the interior of the upper and lower branch of the hyperbola defined 
by equation H3.7|l . Solving equations (|3.8|l for /i and gi and substituting the resulting 



17 



expressions in (|5.1|l we find 



Q{^,t) ^ ^ \ I e*'=^-('=^-^'=)*go(fc)dfc (5.2) 



2n 



^ikx—{k —ik)t 



dD + 



^ik{x- L)-{k^ -ik)t 



dD- 



, , , , N(i-k,t)- N{k,t) 



A(fc) 



A(fc) 

(h Al 

dk, 



dk 



where 



A(fc) = e-^*"^ - e-^(*-'=)'^, (5.3) 

and N{k,t) is defined by (|4.7II . We note that just as with example (c), the determinant 
A(A:) is not zero for k in D. 

Equation H5.2|l shows that there exists a representation of the solution of the Dirichlet 
problem for equation (|1.2|l involving only integrals. It is of course possible to rewrite the 
representation H5.2I) in the series form (|4.12|) using Cauchy's theorem. Indeed, it is easy to 
verify by studying the boundedness of each exponential involved in (|5.2|) that the integrals 
along and dD~ can be deformed to the real line. This deformation passes through the 
hue Im(fc) = i where the zeros km = ^ + 5 of ^(^) li^. Computing exphcitly the residues 
at these poles, and manipulating the result, it is easy to verify that all integral terms cancel 
out. Thus we obtain 

q{^, t) - ^ E e-^'--'"-'" (e^'"" " e'(^-'='")-) (go(fc™) - q,{^ - fc™)) . (5.4) 

mez 

(ii) For the example (a), A(fc) = e**^^ — e~*'^^, thus A(fc) has zeros in D. Hence the part of 
the contour dD'^ along the positive real axis must be deformed to a small circle above each 
of the points km € M., m G Z+. Similarly the part dD~ along the negative real axis must 
be deformed to a small circle below each of the points km S M, to S Z~. 

5.2 The Robin problem for the heat equation 

Let q{x,t) satisfy the heat equation 

qt = Qxx, Q<x < L, t > 0, (5.5) 

with initial condition qo{x) and with Robin boundary conditions at both ends 

q^{0,t) ~ aq{0,t) = hi{t), q,{L,t) - Pq{L,t) ^ h2{t), i > 0, (5.6) 

where a and (3 are given real constants, a ^ (3. 

In this case, X{x,t, k) = qx + ikq, thus the global relation is 

fi+ikfo-e-'''''{gi+ikgo)^qo{k)-e''"*q{t,k), fc e C. (5.7) 

Using the boundary conditions (|5.t)|) . the global relation becomes 

{a + ik)fo~e-'''^{p + ik)go^N{k,t)-e'''*q{t,k), fc e C, (5.8) 
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h2{s) - hi{s)) ds. 



where the known function N{k,t) is given by 

N{k,t)=qo{k)+ f 

Since Ai = — fc, replacing k by — fc in equation H5.8() we find 

(a - i/c)/o - e'^^{(3 - ik)go = N{-k, t) - e^^'q{t, -k), fee C. 
Solving the system (|5.8|) - (|5.1U|) we obtain 

A(fc)/o = e'^^{l3-ik)N{k,t)~e-'^^{l3 + ik)N{-~k,t) 

A(fc)go = {a~ik)N{k,t)-{a + ik)N{~k,t) 



[{a — ik)q{t, k) — {a + ik)q(t, ~k)] 



where 



A(fc) = (a + ifc)(/3 - ik)e'''^ - (a - ifc)(/3 + ik)e-'''^. 
Let fcm be defined by 

km-. (a + «A:™)(/3 - ik^)e'^'"^ = [a - ikrn)W + ikm)e'''''^^ . 



(5.9) 
(5.10) 

(5.11) 
(5.12) 
(5.13) 

(5.14) 



Evaluating equation H5.11|l at A: = km, where A{km) = 0, and using the definition of q{t, km), 
wc find 



_ ifc„)e-*'=-^ - e-*'="^(/3 + ifc™)e*'="'^] (7(x, t)dx 



e-k^t [e^^'"^(/3-iA;™)iV(fc„,t) -e-*'="^(/3 + zA;™)iV(-fc,„,t)] 



(5.15) 



Equation H5.14|l cannot be solved directly for the fc^'s. Thus although the general theory 
implies that equation (|5.15ll can be solved for q{x, t), it does not appear to yield an effective 
representation for the solution of this problem. On the other hand, the general representation 
p.l7|l does provide an effective representation: 



q{x,t) 



^1 


[ e**^ 


x-k^t 


27r \ 


JdD + 





(ik + a)- 



-ikL 



{P + ik)N{-k,t) ~e^''^{l3 ~ik)N{k,t) 



A(fc) 



ahi{t, k) 



dk+ 



^ik{x-L)-k^t 



r,. , ^, {(^ + ik)N{-k,t) - (a - tk)N{k,t) , . 

(ik + — ph2[t, k) 

A(fc) 



dk 



^ikx-k t^^(fc)^;j 



(5.16) 



where dD^ and dD are the boundaries of the domains 

D+ = {\<^rg{k)<^-^}, i?- = < arg(A:) < D = D+ U , 

shown in figure 13 and 



hi{t,k)= / ''hi{s)ds, h2{t,k)^ / e 'h2{s)ds 
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Figure 3: The domain D and the contours dD^ (in red) and dD (in green) for the heat 
equation 



To derive equation H5.16|l we use the solutions H5.11|l . (|5.12|) of the system (|5.8|) - (|5.10ll . In 
the expressions (15.11(1 there appears also a term involving the function q{t,k). This term, 
when multiplied by c*'"'^"'^ *, is analytic and bounded for fc ^ cx) in Z?+. Similarly, the 
term in (|5.12l) involving the function q{t,k), when multiplied by qI-K^-^)-'^ *^ is analytic 
and bounded for fc — )■ cx) in . An application of Jordan's lemma therefore implies that 
these terms give a zero contribution, provided A(fc) ^ for k in D. This is indeed the case, 
and can be seen as follows. The function A(A;) defined by ((5.13(1 is entire and of finite order, 
hence it has infinitely many zeros accumulating at infinity. Moreover, a simple computation 
shows that this function has infinitely many zeros on the real axis, and that all zeros are 
asymptotically on M 10 . It follows that the only one of these zero in Z? is fc = 0, where the 
numerator clearly vanishes. Hence the unknown function q{t^ k) does not contribute to the 
solution representation. 

6 Conclusions 

We have illustrated the applicability of a transform method by solving several concrete 
boundary value problems. It appears that the method is both general and simple to imple- 
ment. Indeed, the only mathematical tools used in this paper are the Fourier transform and 
Cauchy's theorem. 

An effort has been made to minimise technical considerations. In particular, the given initial 
and boundary conditions are assumed to be "sufficiently smooth" . It is possible to work in a 
less restrictive functional class; for problems on the half line x > Q this is done in ^8^, where 
general theorems are proven in appropriate Sobolev spaces. 

The general result about boundary value problems on an interval is that q{x, t) can be 
expressed as an integral in the complex k-plane, see equation (I1.17() . This integral involves 
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the Fourier transform qoi^k) of the initial condition qo{x) and the ^-transforms {fj^Qj} of all 
boundary values {dlq{0, t),dlq{L, t)}. A subset of these boundary values can be prescribed 
as boundary conditions. Thus a subset fp and gq, where p takes N values and q takes n — N 
values, can be computed immediately. The remaining fj's and gj's can be expressed through 
the solution of a system of n algebraic equations obtained from the global relation and from 
the equations derived from the global relation by replacing k with Aj(fc), j = l,...,n — 1. 
The relevant expressions involve the unknown function q{t,k) and the function 1/A(fc), 
where A(fc) is the determinant of the associated system of n algebraic equations. However, 
using the integral representation of q{x,t) (equation H1.17|l l it can be shown that (i): If 
A(fc) ^ in the domain D defined by equation H1.7|l . then the contribution of the terms 
involving q{t, k) vanishes, (ii): If A(fc) has zeros in _D+ (hence also in D~), then the contours 
dD~^ and must be deformed to pass above or below these zeros. For some boundary 

value problems it is convenient to rewrite the contribution of these zeros in terms of a sum 
of residues. In these particular cases there exists an alternative representation consisting 
only of an infinite series. 

From the above discussion it follows that there exists three distinct cases. The basic examples 
Hl.l|l - H1.3|l were chosen in order to illustrate these three cases: 

(a) For equation Hl.l|l . the domain D is the union of the first and third quadrant of the 
complex fc-plane. For Dirichlet (or Neumann) boundary conditions, A(fc) = for k = 

m S Z, thus the contours along the positive and negative real axes must be deformed. 
Furthermore, it can be shown that the relevant integral can be computed entirely in terms 
of a sum of residues. Thus the solution can be rewritten as an infinite series. This is 
consistent with the well known form of the solution as given by a sine series. 

(b) For equation H1.2|l . the domain D is shown in figure For Dirichlet boundary condi- 
tions, A{k) = ioT k = ^ + ^, m G Z, thus A(/c) ^ for A: in D. However, just as in (a) 
above, it is possible to rewrite these integrals in terms of an infinite sum, see equation (|5.2() . 
This is consistent with the solution that can be obtained using the set of eigenfunctions 
corresponding to the complex eigenvalues of the problem, see ^ . 

(c) For equation (jl.8(l . the domain D is shown in figure Q). For the case that the given 
boundary conditions are q{0,t), q{L,t) and qx{L,t), the determinant A(A;) is given by equa- 
tion (|1.20() . In can be shown that the zeros of A(fc) are asymptotically on the three lines 
{k : arg(fc) — 7r/6, arg(A:) = 57r/6 arg(fc) = 37r/2}, and that A(fc) has no zeros in D. 

We note that even in the cases when it is possible to express the integral representation of 
the solution in terms of an infinite sum, the integral form may have some advantages. For 
example, the integral representation, in contrast with the series one, is uniformly convergent 
both at a; = and at x = L. In addition, integrals are more convenient than sums for 
studying the long time asymptotic behaviour of the solution. Furthermore, for many concrete 
examples it is possible to use the integral representation and Cauchy's theorem to compute 
q{x,t) explicitly. It is interesting that in these computations, one does not compute the 
relevant integral by using the residues associated with A(fc) = 0, but one precisely avoids 
these zeros. 

We emphasise that the main advantage of the representation H1.17|l is that the general rep- 
resentation of the solution can be constructed without the need to analyse in detail the zeros 
of A(fc). This is to be contrasted with the classical approach: since the zeros of A(fc) define 
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the discrete spectrum of the associated x-difFerential operator, a detailed characterisation of 
this set is crucial for the derivation of the associated basis of eigenfunctions. This important 
advantage of the method presented here is iUustrated in section 

We conclude with some remarks: 

(1) Coupled boundary conditions can also be analysed using our method. For example, for 
the boundary value problem 

qt + qxxx = 0, q{x,0) = qo{x), g(0,i) = g(L,<) = 0, gj;(L,<) + aga;(0,i) = 0, 

the solution can be expressed in a series of eigenfunctions |12j . This is a consequence of 
the fact that the associated differential operator with the above boundary conditions can 
be made self-adjoint. In this case, it is possible to deform the integrals appearing in our 
representation and to rewrite them in terms of an infinite series. 

(2) Earlier work on two-point boundary value problems for linear evolution PDEs has ap- 
peared in [2], but the existence of zeros of A(A:) was not investigated. The determination 
of conditions for well posedness is presented in [H]. The role of the zeros of the determinant 
A{k) and their determination in some particular cases are studied in ^lOj. 

Appendix 

Let q(x, s) be the Laplace transform of q{x, t), i.e. 

poo 

q{x,s)^ e-'*q{x,t)dt, Re(s)>0. (Al) 

Applying the Laplace transform to equation H1.3|l we find 

qxxx + qx + sq ^ qo{x), Re(s)>0. {A.2) 
The solutions of the homogeneous version of this equation are given by 

q{x, s) = e^"= , A'"^ + A + s = 0. (A.3) 
We distinguish the three roots of the cubic equation (A. 3) by their large s behaviour: 

Xir^-si, X2^-asi, As-^-a^s^, a = e^''*/^. {AA) 
If — 7r/2 < arg(s) < 7r/3 then, for large s, 

-- < arg(s) < -, - < arg(as) < —, — < arg(a s) < —. 
Thus, for large s, 

Re(Ai)<0, Re(Aj) > 0, j = 2,3. (A5) 
A solution of the inhomogeneous equation (A.2) is given by 
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where the constants satisfy the algebraic conditions 

3 3 3 

^/3,=0, 5]A,/3,=0, ^(1 + A2)/3, = 0. (A7) 
11 1 

In the integral 



Jo 



we have x — ^ > 0, thus the integrand is bounded as s ^ oo, if Re(Aj) < 0. On the 
other hand, if we replace with — , then a; — ^ < 0, and the integrand is bounded if 
Re(A_,) > 0. Hence, equations (A. 5) imply that we must choose the following solution of 
equation (A. 2): 

q{x, s) = cie^if^--^) + C2e^^^ + cge^^^ 



where the Cj's are constants independent of x. 

Let /o, go and gi denote the Laplace transforms of the given boundary conditions fo{t), 
ga (t) and gi (t) respectively, 

/•oo fOO POO 

fo{s) = / e-''Mt)dt, ~go{s) = / e-'*go{t)dt, ~gi{s) = / e-''giit)dt. (A9) 
Jo Jo Jo 

The definitions (A.l) and (A. 9) imply 

g(0,s) -/oGs), =go(s), s) - (AlO) 

Using equation (A. 8) to evaluate q{0,s), q{L,s) and qx{L,s) we find the following set of 
three algebraic equations for {cj}^: 

e^^'^ci + C2 + C3 = (3iqo{~Xi) + fo{s), 
ci + e^^^c2 + e^^'^cg = -/32e^^^<7o(-A2) - Pse^'^-qoi-Xz) + go{s), 
Aici + A2e^^^c2 + Aae^^^ca = ^ (32X20^'^ qo{-X2) ~ PiXse^'^^qoi-Xs) + (A.ll) 
where ^ 

qo{X,) = f e^^UoiOdt {A.12) 
Jo 

The determinant of the system (A. 11) is given by 

A = (A3 - A2)e(^^+^-''^i)^ + (Ai - A3)e^^'^ + (A2 - Ai)e^^^. (A.13) 

It can be shown that, as s ^ 00, the zeros of A are on the real axis. Indeed, recall that the 
zeros of the analogous determinant in the complex k plane lie, for large fc, on the half lines 
arg(fc) = 7r/6, arg(fc) = 5tt/6 and arg(fc) = 37r/2. Letting s = i{k^ — k), the determinant 
p.2U|l reduces to the determinant (A. 13), and furthermore the above three rays are mapped 
to the real axis in the complex s plane. 

Since q{x, s) has infinitely many singularities for Re(s) > 0, we cannot define an inverse 
Laplace transform to compute the solution q{x, t) and the method fails. 
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Remarks 



1. Solving equations (A. 11) for the c^'s and replacing the resulting expressions in equation 
(A. 8) we find an expression for q{x, s) which, as s ^ oo, can be shown to be bounded for 
Re(s) > 0, Im(s) 7^ 0. Actually, the relevant proof is identical to the proof of section 0] 
which establishes that for this problem well-posedness demands one boundary condition at 
X = and two boundary conditions at x = L. 

2. Although the classical Laplace transform method fails, it is possible to use the following 
transform instead of the Laplace transform: 



This gives rise to the additional unknown term e^**g(x, t). One must use this term to control 
the singularities of A and then one must eliminate q{x,t) using analyticity arguments. 
3. The application of the Laplace transform involves the solution of two sets of algebraic 
equations, namely equations (A. 7) and (A. 11), while the application of the transform method 
used in this paper uses only the solution of one set of algebraic equations. Furthermore, the 
complex representation is based on the three roots of the cubic equation (A. 3), while the 
representation (|1.17|l is based on the two roots of the quadratic equation 



4. If one uses the transform (A. 14) instead of the Laplace transform, and if one then uses 
the change of variables s — i{k^ — fc), one can then obtain the representation H1.19|l after 
applying certain analyticity arguments to combine terms involving and — to terms 
of the form . 
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